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Abstract 

In this paper we present a dynamical analysis for a coupled tachyonic dark energy with dark 
matter. The tachyonic field <fi is considered in the presence of barothropic fluids (matter and radia¬ 
tion) and the autonomous system due to the evolution equations is studied. The three cosmological 
eras (radiation, matter and dark energy) are described through the critical points, for a generic 
potential V(4>). 
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I. INTRODUCTION 


The role of the tachyon in theoretical physics has a long story. For the most part of 
it, tachyon has been considered an illness whose treatment and cure have been hard to be 
achieved. In the field theory, tachyon is associated with a negative mass-squared particle, 
which means that the potential is expanded around its maximum point. The natural question 
that arises is whether tachyon potential has a good minimum elsewhere. The mexican hat 
potential in the spontaneous symmetry breaking is an example of this perspective. In the 
bosonic string, tachyon appears in its ground state, whereas in supersymmetric string theory 
its appearence can be avoided. However, even in superstrings, tachyon is present in non-BPS 
Dp-branes. However, the potential of such Dp-branes does have a minimum [1, 2] and at 
this minimum the tachyon field behaves like a pressureless gas 0. 

As soon as tachyon condensation in string theory had been proposed, its role in cosmology 
was studied, and plenty of works have thenceforth been done. Concerning the observational 


evidence for the accelerated expansion of the universe 


it, as a dynamical dark energy candidate (see js] for review). Still regarding a dynamical 
dark energy, there exists the possibility of interaction between dark energy and dark matter 

a -—»—™—- 


,S, 


tachyon can be a way to explain 


m 


9, 


When a scalar field is in the presence of a barothropic fluid (with equation of state 
w m = Pm/Pm) the relevant evolution equations can be converted into an autonomous system. 

ihantom 

91 


Such approach was done for uncoupled dark energy (quintessence, tachyon field, p 
field and dilatonic ghost condensate, for instance [6]) and coupled dark energy 


llHl3|. 


Since dynamical systems theory is a good tool to analyze asymptotic states of cosmological 
models, its usage to study coupled tachyonic dark energy is in order. General analyses were 
pointed out in jll. I 2 I] and the case with the potential V(4>) oc <p~ 2 was done in |s|. In this 
paper, we enlarge the previous possibilities considering a generic potential and we assume 
a new form for the interaction between the two components of the dark sector, which in 
turn leads to new fixed points. Then, we solve the dynamical system and we analyze the 
stability of the fixed points, to check what kind of dominated universe can come up. Finally, 
we point out some specific potentials to illustrate the results. 

The rest of the paper is organized in the following manner. In Section [TT] we present 
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the dynamics of the tachyon field, in Section IIIII we analyze the system in the presence of 
barothopic fluids, through the autonomous system. Its fixed points are found and their 
stability analyzed. Section IIVI is reserved for the summary. We use Planck units (h — c — 
1 = Mpi = 1) throughout the text. 


II. TACHYON DYNAMICS 

A. Basics 

Tachyon field 0 is described by the Born-Infeld Lagrangian 


£rt = -^/ = < 7 V r ( 0)\/1 - <^ 0 < 9 M 0 , ( 1 ) 

where V (0) is the tachyon potential. The equation of motion for the tachyon field is 


-V ft d u d> + 


1 + V ( 0 ) 




-dpcpd^cp + 


= 0, 


( 2 ) 


For a homogeneous field, in an expanding universe with Friedmann-Robertson-Walker met¬ 
ric, the equation of motion becomes 


■ V'(<b) 

+ 3 H(f) + = 0. 


(3) 


1-0 2 ' ' ' V(<f>) 

where the prime denotes time derivative with respect to 0. The energy-momentum tensor 
for the tachyon field is 


7 V = -^===4 “ (0) y/l T d a (j)d a (j) : 

Vl + d a <j)d Q (p 

which leads to the energy density and pressure for the fluid given by 


1 — 0 2 

1 - 0 2 . 

With these quantities the equation of state for the dark energy becomes 



W<j> = 


P± 

P<t> 


- 1 , 


(4) 

(5) 

( 6 ) 

(7) 
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thus, the tachyon behavior is between the cosmological constant one = —1) and matter 
one (w<p = 0). Making use of Eqs. (J5]) and (EJ) we get the second Friedmann equation 


“ = _^£L= (1--4? 

“ 3\ji - i 2 ' 2 

From the equation above we see that an accelerated expansion occurs for 0 2 < 2/3. 



( 8 ) 


B. Interacting dark energy 

We now consider that dark energy and dark matter are coupled, but the total energy- 
momentum is still conserved. The continuity equations for both components are 

p<t> + 3i? (p^ + ptf,) — —Q, (9) 


p m + 3Hp m = Q, (10) 

respectively, where Q is the coupling. In principle, the coupling can depend on several 
variables Q = Q(p m , p^, <j), H,t,.. .), so we assume that Q = Qp m p<f,(t)/ H , where Q is a 
constant. With this form, the time dependence of the coupling is implicit in the Hubble 
parameter H. Thus, using (151) . we get an equation of motion for the tachyon similar to ((3j), 
but with an extra term due to the interaction 


+ 3 H cf> + 


V{4>) Q P n 


( 11 ) 


1 -02 1 V(<j>) H 

We are now going to analyze the coupled tachyonic dark energy in the presence of radi 
ation and, of course, matter. 


III. AUTONOMOUS SYSTEM FOR THE MODEL 


The Friedmann equations for the tachyon held in the presence of barothropic fluids are 


H 2 


1 

3 



+ Pm + P 



( 12 ) 
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H 


1 

2 


( vm 2 

\yf^ 


4 

+ Pm + ~P 



( 13 ) 


where the index r stands for radiation. The continuity equations are Eq. CD, for dark 
energy, Eq. (fTOT) for matter and 


p r + 4Hp r = 0 , 


( 14 ) 


for radiation, with H given now by Eq. (TT2l) . To deal with the dynamics of the system, it 
is convenient to define the dimensionless variables [14, Q 


x = 


_ s/Wj 

V s/3H ' 


■JPr 

Z = 


V3H 


A = 


V 1 


r = 


VV" 


( 15 ) 


V 3 / 2 ’ " ” V' 2 ’ 

which are going to characterize a system of differential equations in the form x = f(x, y ,..., t), 
y = g(x, y ,..., t), ..., so that do not depend explicitly on time. For this system 

(called autonomous ), a point (x c , y c , ...) is called fixed or critical point if (/, g,... )\ Xc , yc ,... = 
0 and it is an attractor when (x(t),y(t), ...)—>■ (x c , y c ,...) for t —> oo. 

The dark energy density parameter is written in terms of these new variables as 


o*= p(p 


3 h 2 


( 16 ) 


so that Eq. (Eli can be written as 


1 1(f) + — 1 , ( 17 ) 

where the matter and radiation density parameter are defined by Qi = pi/(3H 2 ), with 
i — m,r. From Eqs. (fT6|) and (fI7|) we have that x and y are restricted in the phase plane 
by the relation 


0 < x 2 + y A < 1, 


( 18 ) 


due to 0 < VLq < 1. In terms of these new variables the equation of state w$ is 


W'f, = x 2 - 1, 


( 19 ) 
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and the total effective equation of state is 


W eff 


Pj+Pr 
P(j> T Pm “I - Pr 


-yV l -1 2 + —, 


( 20 ) 


An accelerated expansion occurs for w e ff < —1/3. 

The equations of motion for the variables x, y, z and A are obtained taking the derivatives 
with respect to N = log a(t), where a(t) is the scale factor and we set the present scale factor 
ao to be one, 


dx 

dN 


(x 


1) 


3x — V3y\ + 3 Q 




dy 

dN 


1 

2 


—VSxyX H— = (x 2 


1 ) + 3 + ^ 2 


dz 

dN 


-2z + | 


3 y 1 


Lv/r 


IX 


X “ 


1) + 3 + ^ 2 


dA 

dN 


-V3Xxy ( T - - 


(21a) 

(21b) 

(21c) 

(21d) 


A. Critical points 


The hxed points of the system are obtained by setting dx/dN = 0, dy/dN = 0, dz/dN 
and dX/dN = 0 in Q2T]) . When T = 3/2, A is constant the potential has the form found 
15] (V(0) oc 4>~ 2 ), known in the literature for both coupled js, 101 and uncoupled 


m 


161 , 


4, 


17] dark energy. The hxed points are shown in Tabled] Notice that y cannot be negative 


and recall that fi r = z 2 . 

The hxed points y c , Xf and yj are shown below 


Vc = 


VA 4 + 36 - A 2 
6 ! 


Q \/Q 2 + 4 


( 22 ) 

( 23 ) 
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(24) 


Vt = 


-A Xf + X 2 x‘j + 12 y/l — x"j 


12(1 


x. 


The point (a) corresponds to a matter-dominated solution, since Q m = 1 and w e ff = 0. 
Points (b) and (c) are radiation-dominated solutions, with f2 r = 1 and w e ff = 1/3. The 
difference between them is that (b) has — 0, while (c) has w$ = —1, which means that 
there is another difference due to their stability, as we will see in the next section. Point 
(d) is a dark-energy-dominated solution with = 1 and w e ff = = —1. Point (e) is also 

a dark-energy-dominated solution (f^ = 1 ) but now the equation of state depends on A, 
which in turn can be either constant or zero. The case with constant A is shown in the Table 
[□and an accelerated expansion occurs for A 2 < 2/\/3. For A = 0 we have y c — 1, x c — 0 
and Wcf, = —1, thus we recover the point (d). All the points shown so far, but (d), were 
also found in [9J, 14, [15j and they do not have any dependence on Q , although the stability 
of the point (e) does, for the interaction showed here. Point (f) is also a matter-dominated 
solution (Q m = 1 and w e ff = 0 ), but the equation of state for dark energy is Wf, = Q 2 — 1, 
leading to an universe with accelerated expansion for Q 2 < 2/3. For this point the coupling 
is restrict to values 0 < Q 2 < 1. 

The last fixed point (g) requires more attention. This solution is valid for Xf ^ Qj and 
for Q / 0, and its behavior depends on Q. In order to have x 2 < 1, we must have Q > 0 
for the case with plus sign in Xf (123]) . while we have Q < 0 for the minus sign case. When 
Q — t|oo| 2 |, Xf — ^0 and yj — > 1 , as it should be, from the restriction (fTSlh Furthermore, the 
fixed point exists for some values of A and Q , due to Eq. (ITSlh In Figured! we show the 
restriciton x 2 + yj as a function of Q, for some values of A. As we see, for A = 0 we have 
yt = 1 only when Xf —> 0 , but this value is never reached since the fixed point (g) is valid 
for Xf 7 ^ 0. Thus, for A < 0 the restriction is not satisfied. From the figure we also notice 
that there is a lower bound for Q, for which the point (g) start existing. In Figure [2] we plot 
yf as a function of Q, for the all the values of A, but zero, shown in Figured!/ 

In Figures |3] and 0] we show the effective equation of state w e ff and the density parameter 
for dark energy O^, respectively, as a function of Q, both for the point (g) (see Table U). 

1 The case for Xf = 0 is the fixed point (d). 

2 Q —i> Too for the plus sign in x/ and Q —> — oo for the minus sign in Xf. 

3 We used the plus sign in Xf (1251) for Figures HI an( I [2] but the minus sign shows similar behavior, with 


Q < 0 due to Eq. m, satisfied for A < 0. The graphics for this case are reflections in the ordinate axis. 
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Notice that depending on the value of Q , the universe does not exhibit an accelerated 
expansion for this fixed point, as we see from the horizontal line in Figured The values of 
Q , at which they reach the upper limits of Eq. (USD and f^, are, of course, the same. This 
feature can also be checked from Figures |T] and HI This fixed point is a tachyonic-dominated 
universe only for a range of values of A with specific Q. We have = 1 for A = 1 and 
Q m 1.35, and for A = 2 and Q ~ 0.35, for instance. The effective equation of state for 
these two values is w e ff = —0.72 and —0.29, respectively. Thus, the former w e ff implies an 
accelerated expansion, while the latter w e ff does not. 


B. Stability around the critical points 

In order to study stability of the fixed points ( x,y,z ) = ( x c ,y c ,z c ), we consider linear 
perturbations Sx, Sy and Sz around them. The perturbations satisfy the following differential 
equations 


d 

dN 


*Sx^ 

5y 


\ Sz J 


/ d£ 

d£ 

JLL 

dx 

dy 

d\z 

dg_ 

dg 

dg 

dx 

dy 

dz 

dh 

dh 

dh 

\dx 

dy 

dz 


5y 
&i 


(25) 


where / = f(x, y,z), g = g(x, y, z) and h = h(x, y, z) are the right-hand side of (121 all . (121bj) 
and (121 cl) , respectively. The 3x3 Jacobian matrix is evaluated at the fixed points and it 



FIG. 1. Restriction xj + yj as a function of the coupling Q, for A = 0,1, 2,4 and 8. The horizontal 
line shows the upper bound for the restriction. Thus, the fixed point does not exist for some values 
of A and Q. 
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A=1 

- A=2 

. A=4 

— - A=8 


FIG. 2. Fixed point yf as a function of the coupling Q , for A = 1, 2,4 and 8. 



A=1 

A=2 

A=4 


Q 


FIG. 3. Effective equation of state for the fixed point (g) as function of Q for A = 1, 2,4 and 8. 
The horizontal line ( w e ff = —1/3) shows the border of acceleration and deceleration. 



FIG. 4. Density parameter of the dark energy as a function of Q, for A = 1, 2 and 4. We show 
in the horizontal line the limit = 1, which does not allow a range of values of Q. 
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posseses three eigenvalues ji \, y 2 and /X 3 . The general solution for the evolution of linear 
perturbations are 


Sx = Ae» lN + Be + Cef* N , 

5y = De^ + Ee + Fe ^ N , (26) 

fc = Ge^ N + He^ + Je M3iV , 

where the capital latin letters (A-I) are integration constants. The stability around the fixed 
points depends on the nature of the eigenvalues, in such a way that they are stable points 
if they have negative values (// 1 , /x 2 , A * 3 < 0 ), unstable points if they have positive values 
(Hi, fj, 2 , fJ >3 > 0 ) and saddle points if at least one eigenvalue has positive (or negative) value, 
while the other ones have oposite sign. We show below the elements of the Jabobian matrix, 
using the notation /) = dj/di 


fx 


2x 


3x — VSyX + 3 Q ( 1 — z 



+ (x 2 



3 Qxy 2 \ 
(1 - x 2 f/ 2 ) 1 


(27a) 


fy = {x 2 - 1) ^-^3A 


6 Qy \ 

v I — X 1 ) 


(27b) 


fz = - 6 Qz{x 2 - 1 ), 


(27c) 


9x 


y 

2 


— V3yX + 3x 


\/l — x 2 


(27d) 


9y 


-V3xy\ - lj/Vl - x 2 + 1(3 + z 2 ), 


(27e) 


Qz = yz, 


(27f) 


3 xz y 2 

= 


(27g) 
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hy = —3yz\/l — x 2 . 


(27h) 


h, = £ 2 


2 + i 


3 y 1 


,(x 2 -\) + 3 + z 2 


, __ . (27i) 

2 |yi — x 2 

Using Eq. (|27|) the eingenvalues of the Jacobian matrix were found for each fixed point 
in Table El The results are shown in Table ED 

he first fixed point (a) is a saddle point, because the eigenvalues //2 and /x 3 have opposite 
signs, while /xi can be either positive or negative. This point represents a matter-dominated 
universe, whose behavior is independent of the coupling Q. Both fixed points (b) and (c) 
lead to a radiation-dominated universe, but the former is a unstable point while the latter 
is a saddle one. The point (d) is a late-time attractor and indicates an accelerated universe. 
The two last eigenvalues of the fixed point (e) are always negative. For this point X 2 y 2 = 3x, 
which in turn is between zero and three, then the first eigenvalue is also negative if Q = 0 or 
Q\ < 0. Therefore, the point describes a dark-energy-dominated universe and can lead to a 
late-time accelerated universe if the requirement /i\ < 0 is statisfied. The effective equation 
of state depends only on A, so the only thing the coupling Q does, in the dynamical system 
analysis, is to change the property of the fixed point. The point (f) is a saddle point with 
a matter-dominated behaviour. However, different from (a), the equation of state for the 
dark energy w<p is no longer zero, but now depends on Q. Finally, the last fixed point (g) 
is stable, since u> e // is between zero and minus one, as it is seen in Table El Thus, /X 2 and 
fi 3 are always negative. The first eigenvalue is also negative because < 1. Therefore, the 
fixed can lead to a late-time accelerated universe, depending on the value of A and Q. Its 
behaviour is shown in Figure 0 for A = 1 and Q = 1.35. 


IV. SUMMARY 

The critical points showed in the last section describe the three phases of the universe, 
namely, the radiation-dominated era, the matter-dominated era, and the present dark- 
energy-dominated universe. The matter-dominated universe can be described by the two 
saddle points (a) and (f), whose difference among them is that the second one has a tachy- 
onic equation of state that depends on Q. There are also two points that in principle can 
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FIG. 5. Phase plane for A = 1 and Q = 1.35. The tachyonic-dominated solution (g) is a stable 
point at x « 0.53 and y « 0.92. 


represent the radiation-dominated era: (b) and (c). The first one is unstable and the second 
one is a saddle point with equation of state for dark energy equals to minus one. A tachyonic- 
dominated universe is described by both points (d) and (e). Point (d) is stable and leads to 
an accelerated universe. Depending on the nature of the coulping and the potential, point 

(e) can also be stable, whose related equation of state depends on A. Finally, the last fixed 
point (g) is stable and can describe an accelerated universe depending on the value of A and 
Q, and it is valid only for interacting dark energy. As we showed in the last section, there 
are options where the last point is viable to describe the present dark-energy-dominated era, 
although its equation of state is different from the observed cosmological constant one 18]. 

Therefore, from the whole analysis, the transition radiation —> matter —> dark energy is 
cosmologically viable considering the following sequence of fixed points: (b) or (c) —> (a) or 

(f) [Q dependent] —> (d), (e) [A dependent] or (g) [Q and A dependent]. 

Although the sequence is viable, the form of the potential dictates whether the fixed points 
are allowed or not. Among several possibilities in the literature, the potential V((p) oc 0 -n , 
for instance, leads to a dynamically changing A (either if A —> 0 for 0 < n < 2, or A—> oo 
for n > 2) 19]. A dynamically changing A is allowed for the fixed points (a)-(d) and (f), 
where for the point (d) A c should be zero. On the other hand, points (e) and (g) require a 
constant A, implying V((p) oc 0 


-2 
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A LLa, llT] Another possibility is a changing A which 
approaches asymptotically to a constant value ( 20 ], thus, the case of constant A may be 
taken into account as an approximation. Overall, there can be a wide variety of potentials 
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where the present dynamical analysis can be applied. 
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TABLE I. Critical points (x, y, z and A) of the Eq. (f2Tj) . The table shows the correspondent 
equation of state for the dark energy ffT9|) . the effective equation of state (1201) and the density 
parameter for dark energy (1161) . for raditation (I2 r = z 2 ) and for matter from Eq. (1171) . 


Point 

X 

y 

z 

A 

Wcj, 

£l(J) 

n r 


W eff 

(a) 

±1 

0 

0 

any value 

0 

0 

0 

1 

0 

(b) 

±1 

0 

±1 

any value 

0 

0 

i 

0 

1 

3 

(c) 

0 

0 

±1 

any value 

-1 

0 

i 

0 

1 

3 

(d) 

0 

l 

0 

0 

-1 

1 

0 

0 

-1 

(e) 

Ay c 

V3 

Vc 

0 

constant or 0 

>?vl i 

3 

1 

0 

0 

W<t> 

(f) 

-Q 

0 

0 

any value 

Q 2 - 1 

0 

0 

1 

0 

(g) 

x f 

Vf 

0 

constant and A > 0 

x}-l 

W eff 

Wj> 

0 1 

W eff 

x f y f \ 

V3 1 





for Q > 0 or A < 0 for Q < 0 







TABLE II. Eigenvalues of the Jacobian matrix in Eq. (1251) and stability of the fixed points. 


Point 

yi 

^2 

M3 

Stability 

(a) 

6(1 ±Q) 

3 

2 

1 

2 

saddle 

(b) 

6 

2 

1 

unstable 

(c) 

-3 

4 

1 

saddle 

(d) 

-3 

-3 

-2 

stable 

(e) 

A 2 y 2 — 3 + V3Q\y 

to 

1 

CO 

AV 2 

2 Z 

stable for Q = 0 or QX < 0 

(f) 

3(Q 2 -l) 

3 

2 

1 

2 

saddle 

(g) 

3 (*> - | 

(xy\ 

Z ) 

3 ( xy\ 4 \ 
2 l C3 3) 

stable 
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